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Abstract 

(N , 

^ . The formalism of flavor conversion of supernova neutrinos is generalized to include 

\ possible differences in the fluxes of the muon and tau neutrinos produced in the star. 

' In this case the radiatively induced difference of the and Vr potentials in matter 

■ becomes important. The and Uj- flux differences can manifest themselves in the 

i effects of the Earth matter on the observed z^e {I'e) signal if: (i) the neutrino mass 

' hierarchy is normal (inverted) ; (ii) the solution of the solar neutrino problem is in the 

LMA region; (iii) the mixing C/e3 is relatively large: jC/esI ^ 10~^. We find that for 
Qh, differences in the ~ Vr {^^ - Vt) average energies and/or integrated luminosities 

^ 20%, the relative deviation of the observed Ve (i^e) energy spectrum at ^ 50 MeV 
^ ' from that in the case of the equal fluxes can reach ~ 20—30% (~ 10—15%) for neutrinos 

crossing the Earth. It could be detected in future if large detectors sensitive to the fe 
(z^e) energy spectrum become available. The study of this effect would allow one to test 
^ ■ the predictions of the z/^, Ur, v^, Vt fluxes from supernova models and therefore give 

. an important insight into the properties of matter at extreme conditions. It should be 

taken into account in the reconstruction of the neutrino mass spectrum and mixing 
matrix from the supernova neutrino observations. We show that even for unequal 
and Vr fluxes, effects of leptonic CP violation can not be studied in the supernova 
neutrino experiments. 
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1 Introduction 



The study of supernova neutrinos is of great interest to both astrophysics and particle 
physics. As far as the astrophysics is concerned, the neutrinos provide precious information 
about the core gravitational collapse, the mechanism of supernova explosion and the prop- 
erties of matter at extreme conditions of large densities and high temperatures. From the 
particle physics point of view, supernova neutrinos allow one to probe the neutrino mixings 
and mass spectrum. Indeed, the properties of the neutrino fluxes are modified by flavor 
conversion effects inside the star (see [|l|]-[|l9l as an incomplete list of relevant works) and in 



the matter of the Earth (see for an early suggestion and |22|, ^ for recent studies). 

Until now, the effects of neutrino oscillations and conversion have been considered in the 
assumption that the fluxes of the muon and tau neutrinos produced in the star are identical. 
Indeed, while the diffusion of Ve in the star differs from that of due to the effect of the 
charged current processes, the non-electron flavors, and Vr-, and their antineutrinos are 
expected to undergo similar interaction effects, dominated by neutral current scatterings. 
Therefore, the non-electron neutrinos have similar fluxes: 

(1) 

In contrast to this, the fluxes of the z/g, Ve and are expected to be rather different, with 
the hierarchy of the average energies {Eg) < (Eg) < {E^) . In calculations of the neutrino 
transport in the star u^, v^-, and v^. are considered equivalent and treated effectively as a 
single species, v^. 

Several effects can break the - and the z/^ - z/^ (or - v^) equivalence, giving rise 
to small differences between the fluxes of these species. Once effects of weak magnetism are 
considered, the and v^, interaction cross section in the matter of the star are different; this 
can lead to a difference between the average energies and luminosities of the z/^ and v^J^ fluxes 
as large as ~ 10% p^]. The production of muons in the core of the neutron star results in 



differences in the transport of z/^ and inside the star, whose effect on the neutrino fluxes 
exiting the star has not been fully explored yet Q 

The possibility of studying these phenomena through the observation of a difference in 
the z/^ and fluxes makes this difference a quantity of particular interest. In addition, 
if an appreciable - difference exists, its effects on the neutrino conversion should be 
taken into account while analyzing the supernova neutrino signals in order to determine the 
neutrino oscillation parameters. 

In this paper we present a generalization of the formalism of neutrino conversion in the 
star and in the Earth to the case of unequal z/^ and fluxes. We address the question 
of whether the z/„ - flux differences could be probed via flavor conversion effects and 



^Another potential source of differences between the v^J_ and fluxes are possible effects of new physics, 
such as lepton number violating processes or violations of lepton universality. Such effects could, in principle, 
influence the whole dynamics of neutrino transport in supernovae, and their detailed study goes beyond the 
scope of the present paper. 
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what the effect of these differences would be on the study of the neutrino mixing and mass 

spectrum with supernova neutrinos. 

The paper is organized as follows. After giving some generalities in sec. H, the effects of 
- Vr flux differences on neutrino conversion in the star are discussed in sec. ^. In sec. ^ 

we study the Earth matter effects and comment on possible signatures of the i/^ - Vr flux 

differences. Discussion and conclusions follow in sec. |[ 



2 Supernova neutrinos and different and z/^ fluxes 



2.1 General features 



At a given time t from the core collapse the original flux of the neutrinos of a given flavor, 
Vai can be described by a "pinched" Fermi-Dirac spectrum, 

L E'^ 

F;^{E,T^,r]^,L^,D) = ^-^^^^^^ ^ , (2) 

where D is the distance to the supernova (typically D ^ 10 kpc for a galactic supernova), 
E is the energy of the neutrinos, is the luminosity of the flavor u^, and Tq represents 
the effective temperature of the z/q gas inside the neutrinosphere. The typical values of the 
average neutrino energies produced by the supernova simulations are 

(Ee) = 10-12 MeV , (Eg) = 12-18 MeV , {E^) = 18-27 MeV , (3) 

and the integrated luminosities of all neutrino species are expected to be approximately 
equal, within a factor of two or so [^6|, The pinching parameter rja takes the values 
?7q, ~ — 2 for and Ur] larger values, r]a ~ — 5, could be realized for the electron flavor 
28| , p7| . The normalization factor F^{ria) is given by 

F3(r/„) = -6Li4(-e'''^) , (4) 

where the polylogarithm function lAn^z) is defined as 

oo 



k=l 



In the absence of pinching, rja = 0, one gets -^3(0) = 77r^/120 ^ 5.68. 
The average energy {Ea) of the spectrum depends on both Tq, and r]a- 



{Ea) = (r{r]a)Ta , (6) 
, , ^Li4(-e''") 
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Eqs. and (J^) give cr(0) = 3.15; larger values of a are obtained for larger (positive) t^q, 
e.g. cr(2) ~ 3.6. 

The luminosity, temperature and pinching of the neutrino flux vary with the time t\ 
if these variations occur over time scales larger than the duration of the burst, r ~ 10 s, 
the integrated flux is still given by the expression (0) with ria — ^ and the total energy 
Ea ^ LaT in place of the luminosity L^,. 

In our study of neutrino conversions inside the star we use the following matter density 
profile of the star: 

P(r) = 10^^ C g■cm-^ (8) 

with C ^ 1 — 15. For p ~ 1 — 10^ g ■ cm~^ expression (|^) provides a good description of 
the matter distribution at any time during the neutrino burst emission and propagation 
through the star |^ |^, For p ^ 1 g ■ cm~^ the exact shape of the profile depends 



on the details of evolution of the star, its chemical composition, rotation, etc. The regions 
at p ^ 10^ g ■ cm~^ are reached by the shock-wave propagation during the neutrino burst 
emission, therefore in these regions the neutrinos cross the shock front, which is characterized 



by large moving density gradients ^1 . 



2.2 and Ur spectral difference 

In the absence of specific predictions for the z/^ - Ur flux differences, which could result from 
a number of physical effects, we consider for these two species the Fermi-Dirac spectra, as in 
eq. (0), with different temperatures, luminosities and pinching parameters. Then the ratio 
of fluxes equals: 

Fr _ (LA (TX F.iv,) exp[E/T, -v,] + l 

FO \lJ \tJ F,{r^^) exp[E/T^ - vA + I ' ^ ^ 

We consider the relative difference 

= - 1 (10) 

and adopt a general parameterization of Ap in terms of the relative differences of the 
temperatures, luminosities and pinching parameters: = (T,. — T^)/T^, = (L^- — L^)/L^ 
and = (77,-— ?7^)/?7^. The presence of charged muons in the star may result in a difference in 
the pinching parameters with no appreciable difference in the luminosities and temperatures, 
while different temperatures or luminosities could appear from the effects of lepton number 
non-conservation or violation of lepton universality. 
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For |er|, |ei|, |e^| <^ 1 the difference Ap can be expanded as 



-4 + 







(1)1 





exp[E/T^-rif,] 
exp[E/Tf, -r]^] + l 
(11) 



for any value of the energy E. It is also useful to consider the high-energy limit, E ^ T^rj^, 
of eqs. d) and (|ll]): 



(12) 



and 



Li3(-e^-) 



Lia 



+ 1 



(13) 



Let us comment on the contributions of e^, cl and to Ap- 

1) . The effect of different luminosities, e/, 7^ 0, gives an energy-independent term in Ap, 
as can be seen in eqs. (P)-(p^). If there are no differences in the temperatures and in the 
pinching factors one merely has = bl. 

2) . If the and temperatures are different, 7^ 0, a critical energy E^'^^ exists at 
which the z/^ and Ut fluxes are equal, Ap = 0. li Ur has a harder spectrum than i/^ (er > 0) 
we have Ai? > at > E^^^ and Ai? < for E < Eq^"^ . For a given e^, the critical energy 
j^corr incrgasgg with the increase of and decreases with the increase of e^. If < the 
opposite situation is realized: Ap > Q E ^ E^^^ and Ap < at E E^'^^. The critical 
energy E^^^ increases with the increase of ep (the decrease in absolute value) and with the 
increase of e^. Above the energy E^?^^ the relative difference A^? of the fluxes increases 
linearly with E in absolute value, as follows from eq. (|13|). Numerically, taking = 7 
MeV, eT = 0.1 and eL = = 0, we flnd Eg'"' ^ 30 MeV and eq. (|Tl|) gives Ap 0.46 
{Ap ~ 0.49) at E = 60 MeV and A^ ~ -0.35 (A^ ~ -0.29) at E = 5 MeV. In general, for 
ep ^ the contribution of the originally non-electron neutrinos to the observed i>e spectrum 
is dominated by the flavor with the softer spectrum in the low energy region, while at high 
energies the flavor with harder spectrum dominates. 

3) . The contribution due to different pinching parameters, 7^ 0, increases with energy 
and turns from negative to positive at ~ (E^) ~ 18 — 27 MeV. At larger energies the 
contribution asymptotically approaches a constant value, as one can see from eqs. (|12D and 
([T3|). This value is small since the ratio Li3(— e^'')/Li4(— e'''") is close to unity and therefore 
it undergoes partial cancellation in eq. ([T3|). Taking 77^ = 1, 77,- = 2 and = = one 
gets Ap ~ 0.14 at E > 40 MeV and A^ ~ -0.46 at E = 5 MeV. 

The description given here applies to the difference of the u^j, and Ur fluxes, F^ and F^, 

can be described in terms of the ratios 



as well. The relative difference Ap 

It 



{Tf - Tfi)/Tp, ep = {Lf - Lfi)/Lfi, and e,, 



= iVf ^ Vfi)/vp. analogously to eq. (|Tl 
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3 Conversion effects in the star 



3.1 Neutrino mass and mixing schemes 

Let us consider a system of three neutrinos, (z/i, z/2, z/3), with masses (mi,m2,m3), related 
to the flavor eigenstates, {u^, v^., Vt)i by the mixing matrix U: Ua = UaiVi- The indices 
a and i refer to the flavor and mass eigenstates respectively. 

We adopt the following parameterization of the mixing matrix (see e.g. 



U = V^H^V^^'V^^ (14) 



where 



= diag{l,l,e'^) , (15) 



and 



C12 S12 \ / Ci3 Si3 \ / 1 

V^^ = I -S12 C12 , V^^=\ 1 , V^^=\ Q C23 S23 

1/ V -Si3 Ci3 / V -S23 C23 



(16) 



with Sij = sin 6ij and Cij = cos 9ij . 

We will consider 3i^-schemes which explain the atmospheric and the solar neutrino data. 
The parameters describing the atmospheric neutrino oscillations are p3| : 

ml-ml = Am^2 = ±(1.5 - 4) • 10~^eV^ tan^ ^23 = 0.48 - 2.1 . (17) 

The two possibilities, Am^g > and Am^g < 0, are referred to as normal and inverted mass 
hierarchies respectively. 

The solar neutrino data are explained either by vacuum oscillations (VO solution), or 
by one of the MSW solutions (LMA, SMA or LOW) with the oscillation parameters 

ml — ml = Am2i, tan^ 612 , (18) 



(see, e.g., numerical values from the analyses [p4 



The angle ^13 is bounded to be small from the results of the CHOOZ and Palo Verde 



experiments [^, |39 

sin^ ^13 ^ 0.02 , (19) 
and the CP-violating phase S is still completely undetermined. 
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3.2 Matter effects: Vra potential; the level crossing pattern 



Due to the very wide range of matter densities inside a supernova, p ~ — 10 g ■ cm" , 
the conditions for three MSW resonances (level crossings) are met. The crossing of and 
Ur levels, the /xr resonance, is driven by the difference of the z/^ and I'r potentials which 
appears at one loop level due to different masses of the muon and tau leptons P0[]: 



Vr 



V ~ 
Vfl — 



r^2 2 



[Up + Tin) In 



w 



(20) 



where Gp is the Fermi constant, m,- and are the masses of the tau lepton and of the 
W boson respectively; Up and Un denote the number densities of protons and neutrons in 
the medium. In eq. (pOf ) the + (— ) sign refers to neutrinos (antineutrinos). The difference 
VrfM is much smaller than the Ve — i^tJ- / effective potential Ve = V^GpUf,, where is the 
electron number density of the medium. One finds 



5 ■ 10-V, 



:2i] 



The conversion induced by the potential ( ^OD is governed by the "atmospheric" mixing 
and mass splitting, Amig and 6*23, eq. (|T^. Depending on the precise value of 6^32 the 
level crossing occurs in the range of densities p ~ — p^r, where p^^ is determined by the 
condition 



Am|2 
2E 



:22) 



Numerically, one gets p^^ ~ 10^ — 10^ g ■ cm~^. In fact, as we will see, the actual position 
of the level crossing is not important. 

The probability P23 of the transition between the eigenstates of the Hamiltonian in 
matter, u^rn and h'2m, can be calculated (following, e.g., |jl8|]) according to 



3X008^923 _ I 



23 



X = -27r 



- ] 

2E 



rfK^(r) 



dr 



(23) 
(24) 



where is the distance from the center of the star to the region defined by eq. 

With the profile (g) and the parameters (|17D, eqs. d^), (||) and (|2|) give x = 500-900, 
corresponding to a very good adiabaticity: P23 = 0. The adiabaticity is good also for steeper 
profiles. We have checked that the adiabaticity condition is satisfied even for the large den- 
sity gradients which could be produced by the shock wave propagation. 



The two other resonances, driven by the z/g 



/ Ur effective potential, Ve, occur at 



much smaller densities, p ^ 10 g ■ cm . For these values of the density the potential Vrf, 
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is negligibly small and it is convenient to consider the evolution of the system in the rotated 
basis of states: 

y'^ = /^V^3t^^ _ (25) 



The second (H-) resonance, 



associated to Am^j 



and 6113, occurs at p = pu 



10^ 



in refs. 11 



10^ g ■ cm~'^. The resonance is in the neutrino (antineutrino) channel if the mass hier- 
archy is normal (inverted). We denote as Ph the probability of the transition between 
the eigenstates of the Hamiltonian in this resonance (the hopping probability). As shown 
Ph increases with the decrease of sin^6'i3, varying from perfectly adia- 
^ 10"'^ to strong adiabaticity breaking {Ph = 1) at 
sin^^ia < 10"^ 

The third level crossing (L resonance) takes place at a lower density, p = ^ 10 g ■ cm~'^. 
It is determined by the "solar" parameters Am^i and 612, eq. ([TB|), and occurs in the neutrino 



batic conversion {Ph = 0) at sin 6^13 -r^ 



channel (see sec. |3.1|) . We denote the hopping probability associated to this resonance as Pl- 
For oscillation parameters in the LMA region the conversion is adiabatic, Pl = 0, while par- 



tial breaking of adiabaticity occurs for other solutions of the solar neutrino problem |11, 22 



The level crossing scheme for the normal mass hierarchy with 6*23 < 7r/4 and a large solar 
mixing angle 612 is given in fig. |l]. The dynamics of neutrino conversions inside the star can 
be considered as occurring in two steps: 

(1) propagation between the neutrinosphere and the layer with an intermediate density 
Pint such that Ph < Pint p^r, ("inner region") and 

(2) propagation from the layer with p ~ pint to the surface of the star ("outer region"). 
Let us consider the evolution in the inner region. At the neutrinosphere the potentials Ve 

and Vrfj, dominate over the kinetic energy differences Am1j2E and the mixing is strongly 
suppressed, i.e. the matter eigenstates coincide with the flavor ones. At the densities 
P ~ Pint, the potential is negligible and the matter eigenstates coincide with the states 
of the rotated basis (|25|) . Using the level crossing scheme of fig. 1 and taking into account 
that the adiabaticity is satisfied in whole range of densities pn < Pint ^ P/^t, we find that 
in the inner region the following transitions take place: 

Z/e->Z/e, I^^^I^r', ^r^V^^, (26) 

provided that the pr level crossing occurs in the antineutrino channel. Therefore in the 
rotated basis {ve-, , ^t') the fluxes of the states (before they enter the H resonance region) 
are 

F{v:) = Fl F{v;) = Fl F{v^) = Fy 

Actually, the observables do not depend on the position of the pr resonance: what matters 
is the flux of the state that crosses the v^, level at the H resonance. 
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Notice that the role of the potential Vr^ is reduced to the suppression of the flavor mixing 
at high densities; the adiabatic change of this potential allows one to relate the fluxes of the 
states of the rotated basis to those of the originally produced flavor states. 

The evolution of the neutrino states in the outer region is determined by the level 
crossings at the H and L resonances. If the hierarchy is inverted, then, as shown in fig. 
^, the fiT and L resonances are in the neutrino channel, while the H resonance is in the 
antineutrino channel. In the inner region the conversions z/^— i>z/^' and v^^v^' take place, 
while the z/g^z/^' conversion occurs in the L resonance. Antineutrinos first undergo the 
v^^v^' and v^—^Vj.' conversions in the inner region. At lower densities, the Ue^i^ti' / i^r 
transitions take place in the H resonance and in the L resonance for large 612 . 

3.3 Conversion probabilities 



As discussed in [|TT1, if = F2 = F2 the effect of the conversion on the flux Fp of electron 



neutrinos observed at the Earth can be expressed in terms of a single permutation parameter, 
p, which represents the survival probability: 

F,=pF° + (l-p)F°. (27) 

We shall generalize now this result to the case F^ ^ F^. 

Let us show that the differences between the t'^ and {p^ and Vr) fluxes can manifest 
themselves only in the following cases: 

(i) normal hierarchy, neutrino channel, 

(ii) inverted hierarchy, antineutrino channel, 

both of them with non-maximal adiabaticity breaking in the H resonance: < 1. Indeed, 
as can be seen from figs. |I] and 0, in the remaining cases (inverted hierarchy, neutrino channel 
and normal hierarchy, antineutrino channel) neutrinos of one of the non-electron flavors are 
completely converted into the third mass eigenstate, for which the electron component 
is small (equal to sin^6'i3). Therefore this flavor does not contribute significantly to the 
observed Uf, (or Up) flux, making the study of the - Vj. inequivalence impossible. The 
same conclusion is obtained in the schemes (i) and (ii) for maximal adiabaticity breaking 
in the H resonance, Ph = 1, since this corresponds to the complete transition of one of the 
non-electron flavors to the third mass eigenstate. 

In what follows we discuss the conversion probabilities for the scenarios (i) and (ii); re- 
sults for the remaining cases can be obtained by a simple generalization. 

Let us first consider the conversion of neutrinos in the scheme with normal hierarchy, 
Am|2 > 0. As we discussed in sec. ^.2|, the neutrinos arrive at the region of the H resonance 
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as incoherent states z/e, f^',i^r' with the fluxes Fj^, Let us consider the further evolu- 

tion of these states. As a result of the adiabatic conversions in the star and/or decoherence 
due to the spread of the wavepackets, the neutrinos leave the star as mass eigenstates. The 
fluxes JF = (Fi, F2, F3) of these states can be expressed in terms of the original flavor fluxes 
F° = F°) according to 



(2^ 



where V is the matrix of the conversion probabilities inside the star in the outer region, 
Via = P{va i^i)- In the approximation of factorized dynamics, (factorization of the 
transition probabilities in H- and L- resonances) the matrix V is fill] 



V 




;i - Ph)Pl 

-Ph)(1-Pl) 
Ph 



(29) 



as can be readily derived from the level crossing scheme, fig. |1|. Then the fluxes of neutrinos 
of definite flavors observed at the Earth, F = [F^ 



F F 



where the matrix S is defined as 



If/. 



are obtained by the projection 

(30) 

(31) 



Combining eqs. (|28| ) -(|31| ) one finds the relation between the original fluxes and the observed 
flavor fluxes: 



F 
Q 



QF^ 
SV . 



(32) 
(33) 



Each entry Q^a of the matrix Q represents the total Va — >■ conversion probability, i.e. 
the probability that a neutrino produced as in the star is observed as in the detector. 
From eqs. 

(0), (H), (HH), dH) and (H), neglecting terms of order sin^^ig, we get 

Q(o) + 2jg« , 



where 



Q 



J = S12C12S23C23S13 cos 5. 



(34) 



(35) 



and 



Q 



(0) 



g(^) = (1 - 2Pl 



Pnh 






(1 - PH)fL 






-Ph) cIJl 


'^23 


;i-p^)(i-^) 


+ sI^Ph 




-Ph) 


■^23 


[l-PH){l-fL) 


+ cI^Ph 













-Ph 1 


-{^-Ph) . 








Ph -1 


{^-Ph) / 









(36) 
(37) 
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In eqs. (|3q) 

fL = PL + il- 2Pl) sin2 9^2 . (38) 
As can be seen from eq. (l37|), the observed Ue flux does not depend on the CP- violating 



phase 5 and, moreover, the terms proportional to J cancel once the sum + Fj- is consid- 
ered. Therefore, since the separate detection of and Ur is impossible, the study of the 
CP-violating phase 6 with supernova neutrinos seems to be not feasible. This result has a 



general validity and is discussed in more detail in sec. 3.4 



We consider now the antineutrino channel in the case of the inverted hierarchy, Am^g < 
(flg. H). Similarly to eqs. (|3^)-(|33[), the fluxes F = {Fg, Ffi, Ff) of the antineutrinos at the 



detector are related to the original fluxes F = (Fg , F?, by 



F = QF , (39) 
Q = SV . (40) 

The matrix V of the i^i conversion probabilities is given by 

/ Ph{1-Pl) Pl {1 - Ph){1 - Pl) \ 
V = PhPl (1 - Pl) (1 - Ph)Pl , (41) 

V {'^-Ph) Ph J 

where Pl represents the hopping probability associated to the L resonance in the an- 
tineutrino channel A form analogous to eq. (p^ is obtained for the Q matrix. 



Q = g(o) + 2JQ(^), where the terms and Q(^) are 

/ PhIl 1-/l {l-PH)fL 

Q(^) = \ d^Pnil - h) + sUl - Ph) c^Jl 4,(1 - Ph){1 ~ h) + sj^Pn 
\ sI,Ph{1 - h) + cUl - Ph) sIJl s^l - Ph){1 - Jl) + cI^Ph 

(42) 

/ \ 

Qii) = _{1_2Pl)\ -Ph 1 -{1-Ph)\. (43) 
V Ph -1 {1-Ph) I 

Here we deflned 

f^ = P^ + [l-2PL)cos^e^2 . (44) 

One can see that the results (p|)- (|43|) can be obtained from those in eqs. (p9D-(^7D by the 
replacement 

^ 1 - Pl , (45) 

or, equivalently, 

fL - II. (46) 
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3.4 Conversion effects and CP violation 



As follows from eqs. (^) and (|43|) , the elements of the second and the third lines of the 
matrices Q^^^ and Q^^^ have opposite signs and therefore the sums + and Fp_ + Ff of 
the observed non-electron fluxes do not depend on S. Let us show that this result holds in 
the general case, without any assumption about the dynamics of neutrino conversion. 

Consider the evolution of the neutrino states at densities p ^ pint, where the potential V^^ 
can be neglected and the conversion can be described in terms of the rotated states (p5D . 
The fluxes of these states coincide with the original Ua fluxes as a result of the adiabatic 
conversion in the pr resonance. From the definition (pSj) it follows that the evolution of the 
rotated states u'^ and therefore the matrix V, eq. (p9l), do not depend on the CP- violating 
phase 6, irrespective of whether the approximation of factorization of transitions in different 
resonances is used. Terms containing 6 may appear in the conversion probabilities only due 
to the projection ( pOD of mass states to the flavor basis by the matrix S, eq. (pT]). We can 
write the sum of the z/^ and Ur fluxes as 

F^ + F^ = J2iS>.i + Sr^)P^aF!^ . (47) 
ia 

As a consequence of unitarity we find that for the 5-dependent parts S^^ of the matrix S 
the following relation holds: 

^2^ + 52^ + ^lf = 0; (48) 

furthermore, in the parameterization (0) we have: 

= . (49) 

The combination of eqs. (^8]) and ( |i9| ) yields 

S^S + S'^ = . (50) 

This means that the sum (|47|) of the and z/^ fluxes is independent of 6. 

Neutrino propagation in an asymmetric matter can give rise to matter-induced T vi- 
olation effects in neutrino oscillations even in the absence of the fundamental CP and T 



violation (i.e. when 6 = 0). In it was pointed out that such effects cannot be observed 
in supernova neutrino experiments as a consequence of the assumption F° = F^ and ex- 
perimental indistinguishability of z/^ and Ur. We note here that this result holds true even 
if the condition Fj^ = F^ is relaxed. This is a consequence of the fact that the effects of 
matter-induced T violation (as well as those of the fundamental T violation) drop out from 
the sum (P7|), irrespective of whether the original z/^ and Uj. fluxes coincide. 
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3.5 Expected spectra 



Let us discuss the observed Ve signal for normal mass hierarchy (the case (i) of sec. ^.31) . 
From eqs. (|32|)-(|37D one finds that the detected flux of electron neutrinos, Fg, is 

Fe = F°g,, + F^QeM + i^°ar , (51) 

= F°Ph/l + F;(1-/l) + F°(1-P//)/l. (52) 

These expressions reduce to the equal fluxes form, eq. (^), if one of the following possibil- 
ities is realized: 

(1) the adiabaticity is maximally violated in the H resonance, Ph = 1 (as discussed in 
sec. |373| ), 

{2)fL = 0, or 
(3) fL = 1. 

As can be seen from eq. (|38|) , the last two possibilities are realized only if the angle 612 is 
small, which is disfavored by the solar neutrino data (see e.g. |5^-|5^). In the following we 
focus on large solar mixing angles and therefore neglect cases (2) and (3). 
Expression (^T|) can be conveniently rewritten as 

= F^PhIl + F'^il - PhIl) + (F° - Fl){l - PH)h , (53) 

where the third term is entirely due to the difference of the z/^ and z/^ fluxes and vanishes if 
either condition (1) or (2) is fulfilled. From eq. ( |55| ) one finds the relative deviation, R, of 
the observed flux from the one in the case of equal fluxes: 

where Ap was given in eqs. (|10|) and (pT]). As follows from eq. (p4D, the sign of R is 
determined by that of Ap. Therefore R vanishes at the critical energy E^^^ and, if e^- > 
{er < 0), we have R < {R > 0) ioi E < E^^"' and R > {R < 0) for E > Fg'''^ 



As discussed in sec. if e-r = e») = no change of sign occurs and the quantity R is 
determined by the difference in the luminosities only. 

(2) 

It is useful to introduce another critical energy F^ defined as the energy at which the 
difference Fg° — F° changes its sign At F <C F^^ the electron neutrino flux dominates 



(2) 

over the z/^ and z/^ fluxes; in contrast, for F ^ F^ the latter fluxes - having larger average 
energy - are dominant. Numerically, taking = 3.5 MeV, = 7 MeV, Le = and 
r]^ = r]f, = one finds E^^^ ~ 20 MeV. 

Taking P^ = (which is the case for solar neutrino parameters in the LMA region), 
from eq. ( p^ we find that 

R 

< — < sin^ 012 , (55) 
Ap 



12 



for any value of the energy E. The upper bound is reached if Ph = and E S> E}j , for which 
the factor in eq. (0) is neghgibly small due to the hierarchy of the temperatures. 

Notice that the condition Ph = corresponds to no contribution of the original z/g flux to 
the observed z/g signal: Qee = (see eq. (p6[)). Thus, the observed Ue flux is given by a 
mixture of the original t'^ and Ur fluxes. The effect is suppressed, i? ^ 0, for Ph ~ 1 or at 
low energies, E <^ Ejj\ where the ratio F^/F^ is large. Numerically, taking P^ = Ph = 0, 
T^ = 7 MeV, er = 0.1, = r/^ = = and tan^ 612 ~ 0.33 one gets i? ^ -0.1 at E = 5 
MeV and R ~ 0.12 at E = 60 MeV. 

In general for ley], I^lI < 0.2 and |e^| ^ 1 the relative deviation R can be as large as 
~ 40 — 50% (in absolute value) at high energies, E Ec\ and ~ 15% in the low energy 
limit, E < E^^\ 

For antineutrino conversion with inverted mass hierarchy the effects of the difference of 
the Ufj_ and Ur fluxes are described by a straightforward generalization of eqs. (^) - (Q). 
The results are similar to those presented here. 



Let us consider the possibility of observation of the difference of the t'^ and u^- fluxes. It 
is easy to see that, despite relatively large values of R found here, the effect of the - Ur 
inequivalence will be very difficult to identify. Indeed, the corresponding distortion of the 
observed Ue spectrum, which for 7^ consists in a broadening of the energy distribution, 
can be mimicked by other effects and therefore does not represent a unique signature of the 
unequal z/^ and fluxes. 

In particular, taking, e.g., eL = = 0, we find that the distorted spectrum is very 
accurately reproduced by the equal fiuxes case with a temperature 

^ (1 - fL)T^ + (1 - PH)fLTr (56) 

of the i^x spectrum. In eq. (^) the z/^ and Ur temperatures are weighted by the same 
coefficients that multiply the corresponding fluxes in the expression of the observed z^e flux, 
eq. (|52|). If > (e^ < 0), we have < < {T^ < < T^) and the differences 
\Tx — T^\ and \T,j. — T^\ are well within the uncertainty of the predicted values of and 
Tr- A numerical example is given in flg. ^ which shows the energy spectrum of events at 
the SNO detector from the z/g +d e~ + p + p detection process with equal and different 
z/^-i^r fluxes. We considered thermal time-integrated spectra as in eq. (|]) with no pinching, 
equal integrated luminosities E^ = E^- = 5 ■ 10^^ ergs and a distance D = 8.5 kpc from 
the supernova. LMA oscillation parameters (tan^6'i2 = 0.33) were assumed, together with 
normal mass hierarchy and Ph = 0. As follows from the flgure, the z/g spectrum obtained 

^Events from the Ve +d e+ + n + n reaction can in principle be distinguished from those from 
Ve +d — > e~ +P+P due to the detection of at least one neutron by capture on deuterium, with characteristic 
time tnd ^ 5 ms, in coincidence with the charged lepton. The coincidence can be established with high 
efficiency (~ 75—80% ^) due to the relatively small event rates of the two processes at SNO (~ 100—200 
events for a supernova at Z? = 10 kpc distance, see e.g. |l^), which guarantees the average time interval 
between two events to be significantly larger than tnd- 
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with different temperatures, = 7 MeV and Tr = 8.4 MeV, is well reproduced by that 
with the common temperature = 7.3 MeV, according to eq. 

Global fits of simulated SNO and SuperKamiokande data Q], performed under the as- 
sumption of equal temperatures of all the non-electron neutrino species and exact equipar- 
tition of the total energy Eb, allow to determine these parameters with precision better 
than ~ 10%. Considering that the data set is dominated by the high statistics in Su- 
perKamiokande, this result should be interpreted as a determination of the effective tem- 
perature and integrated luminosity of the non-electron antineutrinos, and E^. 

It can be checked that, if the and Ut spectra have different pinching factors and 
luminosity difference of ~ 10 — 20%, the equal fluxes case accounts very well for the observed 
Ue spectrum once effective values of and Ex are chosen within their uncertainty intervals, 
as they could be measured by SuperKamiokande. 

In principle, the possibility exists that the differences between the z/^ and Ut (or z/^ and 
u-r) fluxes will be probed at future large detectors (UNO [0, HyperKamiokande etc.) 
once all the neutrino oscillation parameters are known to a sufficient accuracy. 



4 Earth matter effects 

4.1 Unequal fluxes and conversion in the Earth 

If the neutrino burst crosses the Earth before detection, the fluxes of the various neutrino 
flavors in the detector are affected by the regeneration in the matter of the Earth. 

Let us first consider these effects in the scheme with normal mass hierarchy. According 
to a simple generalization of eqs. ( |32D and (|33|), the fluxes F^ = (^F^ , F^ , F:lr') of the 
neutrino flavor states in the detector can be expressed as 



?D 



g®FO , Q® = S®V , (57) 



where the matrix V was given in eq. (p9|). Each entry S®^ of the matrix S® represents the 
probability that a neutrino arriving at the Earth as Vi interacts, after crossing the Earth, 



as a z/q, in the detector: 5® = P(z/j Va)- From eqs. (p9|) and ( pTj) the observed z/g flux 
equals: 

FF = F^Pnff + F;(1 - /®) + F°(l - PH)ff , (58) 

where 

/® = Pl + (1 - 2Pl)S% . (59) 

In eqs. (|5^ ) -(|59| ) terms of order sin^ 6'i3 were neglected 0. 

■^It can be checked that the contributions of these terms to the difference (|o|) are smaller than 10""^, 
and therefore negligible. 
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Combining eqs. (^8]), (|52|) and (|58D -(^9D we obtain the deviation due to the Earth matter 
effects of the observed flux, F^P, from that predicted for conversion in the star only: 

- = (1 - 2PL)freg [PhF', - + (1 - Ph)F^] , (60) 

where the regeneration factor, 

/,,, = 5® -sin2^i2 , (61) 

accounts for the matter effects inside the Earth. It depends on the solar neutrino parameters, 
the neutrino energy, and the direction of the neutrino trajectory inside the Earth. 

In what follows we consider Am^i and 612 in the LMA region. (The Earth matter effect 
is small for the other solutions of the solar neutrino problem; see e.g. the discussion in ref. 
||22|| ). For the LMA oscillation parameters the conversion in the L resonance is adiabatic, 
thus we set = from now on. 

As can be seen from eq. (|60|) , the contribution of F^ to the difference Fj^ — F^ is small 
for large violation of adiabaticity in the H resonance and vanishes for = 1, according to 
the discussion in sec. |3.3| . In this case the difference of the and u^- fluxes clearly has no 
observable effects. Conversely, in the case of good adiabaticity, Ph = 0, the contribution of 
the z/f, flux, F^, is absent and the whole Earth effect is due to the difference F^ — F^. 



It is useful to split the expression ( [60| ) into two terms according to their dependence on 
Ph: 

F^ - = AFj2) + ^F^"^"- , (62) 

where 

AF(2) = U^^Pn (F° - F°) , AFr^ = /,e,(F° - F^) . (63) 

Let us comment on the quantities AFe and AF^""^^. 

(2) 

• AFe describes the whole Earth matter effects if the u^-Ut fluxes are equal. Being 
proportional to the difference F^ ~ it changes the sign at the critical energy F^^ 
discussed in sec. pl5| . If the regeneration factor f^eg is positive (as it is for conversion in 
the mantle of the Earth only |ll], HI) we have AfP > for F < F^^^ and AF]^^ < 



for F > E)j . The quantity AFe is also proportional to Ph and vanishes if Ph = 0, 
according to the fact that if the H resonance is completely adiabatic the flux F^ does 
not undergo the L resonance. Conversely, AFg is maximal if Ph = 1- 



^pcorr jg proportional to the difference of the and z/^ fluxes. As discussed in sec. pT2 
AF^°"' changes sign at F = Fg""*". If freg is positive and > we have AF^°"' < 
for F < E^""^ and AF^°"' > for F > Fg'^''. The opposite signs are realized for 
ex < 0, while there is no change of sign (see sec. [2^ ) if the z/^ and Ur fluxes differ only 
in the luminosity. 
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The relative size of AFe and AF^™'"'' depends on Ph. In particular, in the high energy 
limit, E ^ E^c^ , where F^/Fj^ <^ 1 and the regeneration factor is large ||2^, the condition 
\/\F^orr^^p{2)^ > 1 requires 



Ph<p: 



eq 



peg 



A, 



1 + A, 



(64) 
(65) 



= and 



as follows from eqs. (H) and (|ToD. Taking F = 60 MeV, = 7 MeV 
eT = 0.1 one gets P"^ ~ 0.33. For Ph Z P^ (i-e. Ph ~ 1) the Earth matter effect is 
dominated by the AFp^ ^-c^mr, ,,rV,no \f p., ^ r>^i 



term, while if Ph ^ Ph the larger contribution to the effect is 



given by AF*^' 
fluxes. 

For Ph ^ Ph' 



which, as we have discussed, is entirely due to the differences in the i/^ 



the effect of the u,. 



flux difference is very difficult to identify due 



D 



to the uncertainties in the values of the various parameters: for le-r], I^t] ^ 0.2 the flux F^ 
observed in the detector is very accurately reproduced by the equal v^-Vt fluxes case with 
values of the parameters within their uncertainties. In contrast to this, if Ph ^ P'h the 
term AF^'^'"''' due to the unequal fluxes dominates, and in the limit Ph = it represents the 
only Earth matter effect. The absence of competing effects offers the possibility of probing 
the differences in the - fluxes and therefore makes this case of particular interest. This 
is true especially if Ph is known to be very small (i.e. sin^ 6*13 ^ 10"'^) from independent 
measurements provided, e.g., by neutrino factories or superbeams |^^. A possibility is also 



related to the comparison of charged current and neutral current events from the neutron- 
ization burst at SNO |M. If the value of Ph is unknown, the possibility that the Earth 



effect is due to F^ 7^ F" with Ph = could be distinguished, in principle, from that of equal 
fluxes with a non-zero Ph (eq. (§^)) at least when the two effects have different signs at 
high energies, i.e. for er > and/or ei > 0. 

The description of the Earth matter effects on the conversion of antineutrinos with 
inverted mass hierarchy is analogous to that presented above. Similarly to eqs. (^)-(|6J 
the difference of the observed z/g fluxes with and without Earth effects is given by 



AF 



(2) 



AF^ 



(66) 



AFi') = fr^gPn (F° - F°) , AF-- = /,e,(F° - F°) , (67) 

where 

jl^g = Sg-cos'e,2 , (68) 

and Sfi is the ui conversion probability in the matter of the Earth. The results obtained 
above for neutrinos can be immediately extended to antineutrinos. However, compared to 
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the neutrino channel, the regeneration factor freg is smaller in the antineutrino case, thus 
giving a smaller effect (see e.g. P^). 

It should be noted that, although the potential Vr^ is induced by radiative corrections 
and therefore is much smaller than the potential Ve, it plays a very important role in the 
evolution of neutrino states in the star provided that the original fluxes of u^^ and z/^ are 
different. If one had disregarded it, one would have had an additional suppression factor 
cos 26*23 ^ 1 in the expressions for and AF™^*" and so much smaller Earth matter 

effects. 



4.2 Observable signals 

Let us now discuss the sensitivity of the Earth matter effect to the v^-Vt flux differences. 
Consider the effect in the neutrino channel for normal mass hierarchy, solar neutrino pa- 
rameters in the LMA region (P^ = 0) and sin^ 6*13 ^ 10~^ (i.e. Ph — 0). In this case the 
relative Earth matter effect, r = [F^ — F(.)/Ff,, is entirely due to the differences in the v^-Vr 
fluxes: 

A rpcorr 

r = , (69) 



and, from eqs. (jT^), ( pB| ) and (^Tj) one gets 



A study of the ratio r is presented in figs. in which the cases of different and Vr 

temperatures and luminosities are considered. The effect of different pinching factors, being 
energy-independent at high energy, E ^ i?^-* (see sec. pTT| ), is equivalent to a difference in 
the luminosities and therefore will not be discussed here. 



Let us first consider the situation in which the difference between the and Vr fluxes 
is mainly due to different temperatures and take > and = 0. The ratio r depends 
on the neutrino energy via the regeneration factor f^eg and via the difference A^?, eqs. ([T0| ) 
and (pj]) . The regeneration factor determines the oscillatory behavior of r with the energy, 
while the dependence of r on Ai;' leads to the change of the sign of r at = E(^^^. This 
sign changes from negative at lower energies to positive at higher energies, provided that 
freg > 0. Both freg and Ap iucreasc with energy (see eq. (|11])), therefore the largest effect 
is achieved in the high energy part of the spectrum. It also increases with increasing ct- 
These features are illustrated in fig. ^ which shows r as a function of the energy E for 
different values of er and of the nadir angle 6n. of the neutrino trajectory in the Earth. 
The figure has been obtained using a realistic density profile of the Earth As fig. ^ 
shows, the oscillatory distortions of the neutrino energy spectrum have larger periods for 
larger E, due to the larger oscillation length in matter, and for larger nadir angle 6n (i.e. 
shorter trajectory inside the Earth). The size of the effect increases with the energy due 
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to the linear increase of A^? with E (see eq. ([T3|) ) and to the larger regeneration factor. 
Taking into account that freg ~ 0.7 and that r/ f^eg can be as large as r/freg ~ 0.45 
for E QO — 70 MeV and ~ 0.1, we find r ^ 0.3. For larger temperature difference, 
er ~ 0.2, the effect can reach ~ 70%. 

If < 0, r is positive for E < Eq^^ and negative for E > E^"^^. Besides this change of 
the sign, the general features and the size of the effect are similar to those discussed above 
for positive e-r- This is expected, considering that r is linear in when the latter is small 



see eqs. ([TT|) and ([TOD). The linear approximation is good for le-rl ^ 0.05 — 0.1. 



Let us now consider effect of different luminosities: = and 7^ 0. In this case 
we have Ap = cl (eq. (pT|)), and therefore r ^ ^hfreg, as follows from eq. (^). The only- 
energy dependence comes from the regeneration factor f^eg- For the luminosity differences 
\eL\ ^ 0.2 and f^eg — 0.7 we find |r| ^ 0.15. This is confirmed by the results shown in the 
fig. i 

In general, both the temperatures and the integrated luminosities of the and z/^ fluxes 
may differ. The interplay of these two differences can enhance or suppress the Earth matter 
effect r compared to the cases we have discussed above. Since, as we have shown, the 
larger contribution to the Earth matter effect comes from the temperature difference, it 
is illuminating to consider the effect of a small luminosity difference for a given difference 
e-r of the temperatures. We focus on the high-energy part of the spectrum because at low 
energies the regeneration factor f^eg is small. 

As can be understood from eqs. (H), ([TT|) and (|70[), at high energies the relative difference 
of fluxes Ai?, and therefore the ratio r, is larger (in absolute value) when ez, and er have 
the same sign. 

Conversely, if ei and ex have opposite signs, the effects of the two differences in eq. (jTl|) 
tend to cancel each other. This is illustrated in fig. ^, which shows the ratio r as a function 
of energy for = 0.1 and different values of e^, with the same parameters as in fig. ^. From 
the figures it follows that the enhancement or suppression of the Earth matter effect due to 
the - Ur luminosity difference can be as large as 40 — 50%, depending on the neutrino 
energy. 

If the mass hierarchy is inverted and = 0, the Earth matter effect due to the z/^ - 
iPr flux difference, AF|°''^, appears in the antineutrino channel without competing effects 
{afP = 0, see sec. [4.1|) . The expression for the relative Earth matter effect r can be 
obtained from eqs. (^), (0) and (|66D -(|67|): 

-freg • (71) 



As can be seen from eqs. (^) and (ffTD, for equal relative differences of fluxes, Ap = Ap, 
the quantity f differs from r by small terms of order A|n and by the suppression factor 
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fregl freg < 1- Foi the antineutrino channel plots analogous to those in figs. | and ^ are 
presented in figs. and ^. 



As a further illustration, in figs. and 10 we show the binned energy spectra of the 



events expected in the heavy water tank of SNO for different and v-r fiuxes with and 
without Earth crossing. In the calculation of the expected numbers of events the detection 
efficiency and energy resolution have been taken into account following the discussion in ref. 
22| . In fig. ^we considered unpinched spectra with = 0, e-p = 0.1 and the other neutrino 



parameters the same as in fig. ^. The nadir angle On = 0° was taken; the distance D and 
the integrated luminosities are the same as in fig. ^. Fig. [l^ shows the same as fig. ^ but 
with Arn^^ = 3 ■ 10"^ eV^ and 6^ = 40°. 

As follows from the figures, the Earth matter effect caused by the difference of and 
Ur fluxes leads to a distortion of the observed energy spectrum as a whole or in isolated 
energy bins, depending on the depth of the neutrino trajectory in the Earth. This distor- 
tion amounts to 10 — 20% at most, and, due to the expected limited number of events, is 
not statistically signiflcant. This deviation could have the signiflcance of 3a or larger for 
a supernova at a smaller distance (2-3 kpc) and/or in a larger experiment. For instance 
the proposed liquid argon detector LANNDD is expected to observe about 3000 events 
49| - more than 10 times larger statistics than SNO - for the parameters used in flgs. ^ 



and |1^. The sensitivity of this detector to the Ue energy spectrum, however, remains to 
be investigated. Moreover, as discussed in |2^, establishing the Earth matter effects may 



require the comparison of the spectra observed by the detectors at different locations, and at 
present it is not clear if more than one large z/e detector with the necessary characteristics^ 
will be built in future. 



The Earth matter effect in the antineutrino channel, eq. (|71|) , can be detected via the 
t'e +P ^ + n reaction in Cherenkov or scintillator detectors The distortions of the 
observed energy spectrum compared to the case of no Earth crossing by neutrinos do not 
exceed 10%. Establishing this difference with at least 3a statistical signiflcance would re- 
quire a comparison of data from two large detectors. We flnd that, unless the distance 
to the supernova is very small {D =2 - 3 kpc), two detectors of SuperKamiokande (SK) 
size do not have enough sensitivity. In contrast to this, large water Cherenkov detectors of 

301 



next generation like UNO g5|, HyperKamiokande pl or TITAND g (~ 18, 40 and 60 
times larger than SK respectively) could establish the effect with signiflcance 3a or larger 
for distances up to ~ 10 kpc at least. 



The results for the time-integrated energy spectra presented here can be extended to 

^Besides the high statistics, the observation of the Earth matter effect requires (i) separate detec- 
tion of Ve and z/g + t'/i + i^t signals, (ii) separate detection of neutrinos and antincutrinos, and (iii) the 
reconstruction of the energy spectrum. 

^We checked that the contributions of charged current scatterings on nuclei, e.g. i^e +^^0 -^^^ F + e^, 
to the event rates is smaller than 2 — 4% and therefore negligible. 
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the spectra of events collected in any (small) time interval within the duration of the burst. 
Taking into account that the difference between the fluxes of the original non-electron neu- 
trinos may vary with time, a study of the time dependence of the signal would be of interest 
to test the effect we have discussed. This is especially true if a specific time dependence of 
the z/^ - Ur flux difference is predicted by supernova models. 



5 Discussion and conclusions 

If the solution of the solar neutrino problem is in the LMA region, the Earth matter effects 
appear in the spectra of v^. or from a supernova or in both, depending on the type of the 
neutrino mass hierarchy, mixing angle ^13 and possible differences between the t'^, Vt-, and 
Vj. fluxes. In particular, if (i) = (F? = (ii) the hierarchy is normal (inverted), 
and (iii) the conversion in the H resonance is adiabatic, Ph — (i.e. sin^ ^13 > 10"^), the 
Earth matter effects exist in the Ue {I'e) channel only |Tl|, |22l . 



This is no longer true if there are differences between the original z/^ and Ur or and 
9r fluxes. If F° 7^ F^ {F^ 7^ F^) we find that, even for = 0, effects of regeneration in 
the matter of the Earth as large as ~ 70% (~ 20%) can take place in the {ve) spectrum 
for normal (inverted) hierarchy. 

Therefore, if Ph is known to be very small from independent measurements (e.g. from 
neutrino factories or superbeams), the observation of a large Earth effect in the antineutrino 
channel and of a smaller one in the neutrino channel will testify for normal mass hierarchy 
and differences between and Vr fluxes. Moreover, the character of the - difference 
can be found from the sign of the observed Earth effect: a positive (negative) effect at high 
energies will correspond to > or = with > {T^ < or = with 
Lr < L^). 

If the same experimental result is observed but Ph is unknown, an ambiguity exists 
between the two possibilities: (1) F° ^ F° with Ph = 0, and (2) F° = F° with < < 1. 
In some circumstances the ambiguity can be resolved by considering the sign of the effect: if 
the effect is positive the case (1) will be established, with > or = with > L^. 
Moreover, the information that Ph <C P^ is obtained. If the effect is negative, the case (1) 
is allowed only if < or = with < L^. Therefore a full discrimination requires 
the knowledge of the sign of and e^. 

If a large Earth matter effect is observed in the spectrum together with a small one 
in the Ue spectrum, the inverted hierarchy and the existence of the t'^ - Ur flux differences 
will be established, with implications analogous to those given here for the neutrino channel. 

Let us summarize our main results. 

1. The formalism of flavor conversion of supernova neutrinos has been generalized to the 
case of unequal fluxes of the non-electron neutrinos. A proof was given that, even for dif- 
ferent fluxes of z/^ and neutrinos, the effects of CP violation are not observable in the 
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supernova neutrino signal since they do not affect the observed z/g ffux and cancel in the 
sum of the z/^ and Vr fluxes in the detector. 

2. It was found that possible - {p^ - v^) flux differences modify the observed Vf, {p^ 
flux compared to what is predicted for equal fluxes if the mass hierarchy is normal (inverted) 
and the adiabaticity violation in the H resonance is not maximal, Pu < 1. 

3. For conversion in the star only, the effect of these differences consists in a distortion 
(broadening) of the observed z/g or z/g energy spectrum. For relative differences between the 
temperatures and/or integrated luminosities of the muon and tau neutrino fluxes ^ 20% the 
relative deviation R of the spectrum can be as large as i? ~ 0.4 at high energies. However, 
the spectrum is consistent with the undistorted one within the uncertainties due to the 
approximate knowledge of the neutrino temperatures and luminosities. 

4. If the solution of the solar neutrino problem is in the LMA region and the neutrino 
trajectory crosses the Earth, oscillatory distortions of the u^. or Uf, energy spectra occur. 
For Ph ~ 1 ( sin^ < 10~^) the effect of the z/^ - u^- {u^ - z/,-) flux differences is small 
and indistinguishable from the equal fluxes case due to the uncertainties in the neutrino 
temperatures and luminosities. 

5. Conversely, if Ph — (sin^ ^13 ^ 10~^) and the hierarchy is normal (inverted), the Earth 
matter effect in the neutrino (antineutrino) channel is entirely due to the - [p^ - u^-) 
flux difference and therefore the existence of the flux difference can, in principle, be estab- 
lished. In the neutrino channel, the relative deviation r of the observed spectrum from the 
prediction in the case of no Earth crossing can be typically as large as r ^ 0.2 — 0.3 at high 
energy, ^ 70 — 80 MeV, and can reach r ~ 0.7 for the largest considered values of the 
difference Tr — T^. This effect could be observable in future if large Ue detectors with good 
energy resolution become available. In the antineutrino channel the effect does not exceed 
~ 10 — 15% and could be tested by the next generation of large Cherenkov detectors. 

6. The observation of the Earth matter effect that we considered here would be a signal of 
the difference between h>^ and z/^ (or z/^ and Ur) fluxes, with interesting implications for the 
models of neutrino transport in the star and for the underlying physics. Jointly with other 
observations, it would give information on the neutrino mass hierarchy and on the mixing 
angle 9is. The effects of the possible - u^. flux difference should be properly taken into 
account in the interpretations of the results of the supernova neutrino experiments. 

7. The potential Vr^ is induced by radiative corrections and therefore is much smaller than 
the potential Ve. Nevertheless, it plays an important role in the evolution of the neutrino 
states in the star provided that the original fluxes of z^^ and Vj. are different. Had one 
disregarded it, one would have arrived at a much smaller Earth matter effect, suppressed 
by the factor cos 26*23 1- 
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Figure 1: The level crossing scheme of neutrino conversion in the matter of the star. The 
figure refers to the case of normal mass hierarchy with 623 < tt/A and a large solar angle 612- 
The semi-plane with negative density, Ug < 0, describes the conversion in the antineutrino 
channel. 
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Figure 3: The energy spectrum of the u^. +d e~ + p + p events expected at SNO for 
conversion in the star only with equal (dotted and dashed curves) and different and u^- 
fluxes (solid curve). The dotted curve was obtained considering normal mass hierarchy with 
Ph = 0, = 7 MeV, r]^ = 0, and LMA oscillation parameters (P^ = and sin^ 26*12 = 0.75). 
The solid curve refers to the same parameters with the difference = 0.2. The dashed 
curve mimicks the effect of the different v^-Vt fluxes with = Tr = 7.3 MeV. A distance 
D = 8.5 kpc from the supernova and integrated luminosity = E^- = b ■ 10^^ ergs were 
taken. 
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Figure 4: The relative Earth matter effect in the neutrino channel, r, as a function of the 
neutrino energy for different values of the trajectory nadir angle, and for a number 
of values of the relative - v^- temperature difference ex- We considered normal mass 
hierarchy, Ph = 0, = 0, = 7 MeV, sin^ 2^12 = 0.75 and Am^^ = 5 ■ lO'^ eV^. 
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Figure 5: The same as in fig. ^ but for = and different values of the relative differences 
of luminosities, ei. 
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Figure 6: The relative Earth matter effect in the neutrino channel, r, as a function of the 
neutrino energy for = 0.1 and different values of e^. We have taken 9^ = 0°; the values 
of the other parameters are the same as in fig. ^. 
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Figure 7: The same as in fig. H but for tlie antineutrino cliannel and inverted mass liierarcliy. 
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Figure 8: The same as in fig. |^ but for = and different values of the relative differences 
of the integrated luminosities e^. 
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Figure 9: The spectrum of events expected at SNO from the Ue +d — > e~ + p + p process 
with and without Earth crossing. We took = and = 0.1, a distance to the supernova 
D = 8.5 kpc, nadir angle 9n = 0° and integrated luminosities E^j^ = Er = 5 ■ 10^^ ergs. The 
other parameters are the same as in fig. ^. 
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Figure 10: The same as fig. || but for Am^i = 3-10 ^ eV^ and 9n = 
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